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ABSTRACT 



Extensive grids of photoionization models have been 
calculated for single star H II regions and evolving 
starbursts. We illustrate the predictions for IR fine 
structure lines which are used to analyse the stellar 
content, and derive properties such as the age and 
IMF. The impact of recent ionizing fluxes on the IR 
lines are shown. First comparisons of our starburst 
models with IR-diagnostics and the ISO observations 
of Genzel et al. (1998) are also presented. 
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1. INTRODUCTION 



The analysis of fine structure lines is crucial to inter- 
pret IR observations of ultra-compact H n regions, 
giant H II regions and starburst galaxies in terms 
of their stellar content. Through photoionization 
modeling the nature of the ionizing sources of H n 
regions and properties of starburst regions such as 
the age and the IMF (slope, upper mass-cutoff) can 
in principle be determined. Considerable progress 
has been made in the recent year in this respect 
through the use of new stellar fluxes, spectral energy 
distributions appropriate for evoling stellar popula- 
tions at different mctallicitics and including detailed 
observational constraints (e.g. Garcia- Vargas et al. 
1995, 1997, Stasinska & Leitherer 1996, Stasinska & 
Schaerer 1997, Gonzalez-Delgado et al. 1998). None 
of these studies have, however, focusscd on IR ob- 
servation which are now becoming widely available 
thanks to ISO. 

In this contribution we present exploratory results 
regarding predictions of IR lines from two sets of 
recent photoionization models: 1) Single star H II 
regions using the latest ionizing fluxes from non- 
LTE, line blanketed atmospheres including stellar 
winds (CoStar models, Schaerer & de Koter 1997) 
and, 2) Photoionization models for evolving star- 
burst populations based on the Leitherer & Heck- 
man (1995) synthesis models. Detailed results from 
new photoionization models for starbursts including 
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the CoStar models and updated stellar tracks will be 
presented elsewhere. 



2. SINGLE STAR HII REGION MODELS 



In Figure 1 we illustrate the behaviour of the line 
ratios of some important fine structure lines observ- 
able in the SWS range and with ISOCAM/CVF. This 
figure, based on the model calculations of Stasinska 
& Schaerer (1997, hereafter SS97), shows the depen- 
dence of the line ratios on the stellar temperature 
and - importantly - also on the ionization parame- 
ter. The comparison between the top and bottom row 
shows the differences obtained using either the recent 
CoStar atmosphere models or Kurucz plane parallel 
LTE models. Due to the harder spectrum predicted 
by the CoStar models above ~ 35-40 eV (Schaerer & 
de Koter) a higher excitation is generally predicted 
for a given T e g than with Kurucz models (see SS97). 

Evidence supporting such an increase of the hardness 
of the ionizing flux comes from the KAO observations 
of Rubin and coworkers (cf. Simpson et al. 1995). In- 
deed, the CoStar models including stellar winds, non- 
LTE effects and line blanketing are thereby able to 
provide a satisfactory description of the Ne ++ /0 ++ 
ionic fractions (see SS97, also Sellmaier et al. 1996). 

Encouraged by these successes atmosphere models 
including the most relevant physical processes (stel- 
lar winds, non-LTE, blanketing) continue to be im- 
proved. See e.g. Pauldrach (1997) and Schaerer 
(1997) for recent progress reports. In addition more 
stringent tests using ISO observations of nebulae with 
well known ionizing sources would be extremely valu- 
able. Detailed photoionization modeling including 
constraints on the nebular geometry and likely also 
their dust content will be necessary to achieve such 
a goal. 

In Fig. 2 we plot the predictions for single star H n 
regions for some quantities of interest for the inter- 
pretation of integrated populations discussed later. 
The middle and right panels show the ratio of lumi- 
nosity in the Lyman continuum (assuming an aver- 
age photon energy of 16 eV) L^ yc to the [Ne n] and 
[S in] line luminosity respectively. It is useful to un- 
derstand the behaviour of the [Ne n]/[S ill] 12.8/33.5 
ratio (left) since either of these lines is sometimes 
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Figure 1. Top row: Photoionization models using CoStar fluxes (dwarf models for T e e> 35000 K, supergiant 
models otherwise. Bottom: same with Kurucz models. Models at solar metallicity with varying ionization 
parameter U over 4 dex (typically \ogU ~ -5 to -1). Solid lines connect models with U typical for giant H II 
regions or for bright single star H II regions ("reference model" of SS97). 



used as a proxy in case the other line is unobserv- 
able (cf. Genzel et al. 1998). Indeed for a given T e g 
one obtains a relatively small variation (~ 1 dex) of 
[Ne n]/[S in] with the ionization parameter. In all 
our "reference models" (constructed assuming a gas 
density of 10 cm -3 ) we predict [Ne n]/[S in] < 1. For 
stars with T e g > 36 kK the differences between the 
CoStar and Kurucz fluxes become non-negligible for 
[Ne n]/[S in]. Some implications of these finding are 
briefly discussed below. 

3. PHOTOIONIZATION MODELS FOR 
EVOLVING STARBURSTS 



IR fine structure lines provide a new tool to study the 
hot star population of obscured starbursts. Through 
comparisons with photoionization modeling this al- 
lows in particular to constrain the star formation his- 
tory and the upper end of the initial mass function 
(e.g. Lutz et al. 1998). Empirical diagnostic diagrams 
based on IR lines also allow the distinction between 
starburst and AGN powered sources (Genzel et al. 
1998). Theoretically these diagnostics remain, how- 
ever, to be understood and quantified. 

To address such questions photoionization models 
representing evolving starbursts at different metal- 
licities are nowadays required. In the following we 
will briefly illustrate some predictions from such 
models for fine structure lines in the ISO LWS 
range. The predictions shown here are drawn 
from the models of Stasinska & Leitherer (1996, 



hereafter SL96, available by anonymous ftp from 
ftp.obspm.fr / pub / obs /grazyna/cd-crete) . Detailed 
results from new photoionization models for young 
starbursts based on the most recent evolutionary syn- 
thesis models of Schaerer & Vacca (1998) including 
the latest Geneva stellar tracks and CoStar model 
atmospheres (cf. above) will be presented elsewhere. 

In Figure 3 we show the computed temporal evolu- 
tion of the line ratios of [Ne m]/[Ne n] 12.8/15.5 
/zm, [Ar m]/[Ar n] 8.9/6.9 fim, and [S iv]/[S ill] 
10.5/18.7 /im (all top row) for an instantaneous burst 
at 1/4 solar metallicity with a Salpeter IMF and 
M up = 100 M surrounded by a gas cloud of den- 
sity n=10 cm" 3 . Different symbols represent models 
with different (initial) ionization parameters: typi- 
cally one has \ogU ~ — 5 (circles) to logt/ ~ — 1 
(crosses). Note that in photoionization models of 
evolving starbursts, the dependence of U with time 
is a function of the adopted geometry, as stressed by 
Stasinska (1998). In the models of SL96, U decreases 
slightly with time. 

As expected the excitation decreases with the age of 
the burst. It must, however, be reminded that the 
considered lines also strongly depend on the ioniza- 
tion parameter (see Fig. 3). Indeed the line ratios 
shown here vary typically by three orders of magni- 
tudes for variations of U of 4 orders of magnitudes. 
For comparison a change of the upper mass-cutoff 
from 50 to 100 M Q increases the same line ratios by 
less than one order of magnitude (e.g. Lutz et al. 
1998). In general U is poorly constrained. Although 
optical spectra of H n galaxies can be reproduced 
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with models covering only a fairly small range of U 
(~ 1 dex, SL96) this is not necessarily the case for 
ISO observations of starburst galaxies. This question 
needs a careful analysis before meaningful results can 
be derived on the IMF, age etc. 

Genzel et al. find an empirical value of 
[Ne n] /[S ill] = 1.7 ± 1.3 (or 1.2 ± 0.6 excluding 
the two largest outliers) from their sample of 10 star- 
bursts. The SL96 models constructed with a density 
of n=10 cm -3 typically predict [Ne n] < [S in] (this 
is true also for solar metallicity models, not shown 
here). However, the SL96 models with n= 10 4 cm~ 3 
reach [Ne n] / [S in] of 1.5, due to collisional deexci- 
tation of [S in]. The electron densities deduced from 
[S in] 18.7/33.5 for a subset of starburst galaxies from 
the Genzel et al. sample are typically < 10 2 ~ 25 , so it 
seems that collisional deexcitation is not sufficient to 
explain the large [Ne n]/[S in] observed. Among pos- 
sible factors to reduce the discrepancy are an Ne over- 
abundance (which however needs to be understood in 
more detail), or a softer ionizing radiation field. We 
conclude that the empirical value [Ne n] / [S in] ~ 1.7 
for starbursts used in the diagnostic diagrams of Gen- 
zel et al. (1998) is theoretically not well understood 
yet. 

For comparisons of the Lyman continuum luminos- 
ity Lhyc to the bolometric luminosity of starbursts, 
ULIRGS, etc. it is interesting to be able to derive 
Lhyc from the IR fine structure lines. This requires, 
however, some empirical calibration or the use of 
models. "Empirical" values of L-^ yc / L([Ne n]) = 
64 ± 37 and L Lyc /L([S in]) = 105 ± 61 are obtained 
from the starburst sample of Genzel et al. excluding 
NGC 5253, which shows much larger values. It must 
be reminded that L^ yc is derived from the Brackett 
recombination lines assuming case B and an average 
photon energy of 16 eV, and neglecting the absorp- 
tion of Lyc photons by dust. 



4. CONCLUSIONS 



We have illustrated the impact of the recent CoStar 
ionizing fluxes on IR fine structure lines for single 
star H II regions and shown predictions of IR lines 
for evolving starbursts resulting from calculations of 
extensive grids of photoionization models (see SL96, 
SS97). 

As shown earlier, ionizing fluxes from non-LTE, line 
blanketed atmospheres including stellar winds pre- 
dict a harder spectrum (Schaerer & de Koter 1997) 
leading to higher excitation in the nebula. These pre- 
dictions are supported by comparisons with optical 
and pre-ISO (KAO) IR observations of H II regions 
(SS97 and references therein). Detailed models of 
ultra-compact H n regions or H II regions with well- 
known stellar content based on ISO observations and 
other observational constraints will be extremely use- 
ful to better constrain the atmosphere models and 
probe the importance of other effects (e.g. dust) usu- 
ally neglected in photoionization models. Such stud- 
ies also provide important tests for our understanding 
of more complex objects like starbursts. 

For starbursts we have illustrated the behaviour of 
line ratios which can in principle be used to constrain 
burst properties such as the age, IMF slope and up- 
per mass-cutoff. We have also studied the behaviour 
of the [Ne n] / [S m] ratio and diagnostics of L^ yc from 
fine structure lines extensively used by Genzel et al. 
(1998). We have shown that a thorough analysis of 
the properties of the nebular gas (ionization param- 
eter, geometry, dust content, total amount of gas) is 
needed for reliable determinations of the stellar con- 
tent of the nebulae observed by ISO. 
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Figure 2. Same conventions as Fig. 1. 
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Figure 3. Photoionization models for evolving starburst (instantaneous burst, Sal-peter IMF, M up —100 M Q , 
Z = 1/AZq) and varying ionization parameters (log U ~ -5 to -1, circles to crosses). 



